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ABSTRACT 


A  two-dimensional  planer  labyrinth  seal  test  rig  was  designed  to  operate  with  air 
supplied  at  45  psig  and  temperatures  up  to  150  °F.  The  rig  q;)erates  with  a  manually 
specified  test  section  pressure  up  to  30  psig  yielding  Mach  numbers  to  0.9  and  gap 
Reynolds  numbers  to  100,000.  The  air  flow  rate  through  the  seal  vdll  be  controlled  by 
setting  inlet  pressure  and  adjusting  an  oudet  control  valve.  The  test  section  measurements 
are  18  inches  wide  by  1.5  inches  depth  by  6  inches  in  length  and  provides  for  10:1  large 
seal  "  geometry  seals  to  be  used  to  facilitate  measurements.  Design  maximum  seal  gq>  size 
is  0.15  inches.  The  test  section  has  a  glass  viewing  port  to  allow  flow  field  measurement 
by  non-intrusive  means  such  as  Laser  Dopplo-  Velocimeter  (LDV)  with  seals  containing  up 
to  5  sealing  knives.  Measurements  of  pressure,  temperature  and  flow  fields  can  also  be 
simultaneously  measured  by  probes  inserted  in  the  seal  itself,  or  mounted  on  the 
removable/replaceable  top  plate.  Inlet  flow  is  conditioned  through  the  use  of  a  dump 
diffuser  incorporating  screens,  honeycombs,  expansion  and  contraction  portions.  The  inlet 
flow  to  the  test  section  can  be  modified  from  uniform  to  various  non-unifearm  conditions  by 
employing  proflle  generators  such  as  screens  and  winglets.  A  detailed  mechanical  design 
has  been  conducted  including  stress  analysis  and  seal  flow  rate  predictions. 
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I.  INTRODUCTION 


There  has  been  rapid  development  of  increasingly  more  efficient  and  power  dense 
turbomachinery,  particularly  in  the  area  of  gas  turbines.  The  U.S.  Navy  is  particularly 
intei'ested  in  improving  turbomachinery  since  most  of  its  active  surface  combatants  use  gas 
turbines  for  main  propulsion  and  electrical  power  generation.  In  addition,  every  aircraft  in 
the  Navy  inventory  is  driven  by  some  form  of  gas  turbine.  Higher  efficiency  power  plants 
directly  translate  to  better  range  and  longer  on-station  time  without  refueling,  and  lower  life 
cycle  operating  costs.  Higher  power  density  generally  results  in  smaller  machinery  and 
increased  performance  (speed  and  acceleration) . 

Most  of  the  recent  advancements  in  gas  turbine  efficiency  have  been  realized 
through  increased  thermodynamic  efficiencies  of  the  components,  higher  cycle  pressure 
ratios,  and  higher  turbine  inlet  temperatures.  Signif  cant  resources  are  continuing  to  be 
directed  toward  improved  turbomachinery.  For  example,  a  jet  engine  manufacturer  may 
spend  tens  of  millions  of  dollars  and  years  of  research  in  an  attempt  to  improve  the 
efficiency  of  a  gas  turbine  compressor  by  only  a  fraction  of  a  percent.  In  most  cases  the 
risk  level  is  very  high  since  there  is  no  guarantee  that  further  development  will  improve  the 
efficiency.  While  most  of  the  losses  in  efficiency  are  due  to  the  losses  generated  in  the 
primary  air  flow  path,  some  interstage  losses  occur  as  a  result  of  ineffective  sealing 
between  higher  and  lower  pressure  sides  across  each  stage.  Non-contact  labyrinth  seals  are 
typically  employed,  and  are  one  method  to  minimize  the  interstage  losses  (and  improve 
efficiency)  within  a  gas  turbine  by  reducing  this  high  to  low  pressure  leakage.  Figure  1.1  is 
a  two  dimensional  representation  of  a  3  dimensional  annular  labyrinth  seal.  The  object  of 
this  seal  is  to  minimize  flow  rate  by  successively  throttling  the  flow  in  one  or  more 
chambers  to  destroy  kinetic  energy.  Most  approaches  to  increasing  seal  effectiveness 
pertain  to  closing  the  sealing  gap,  adding  more  "torturous"  passages,  and  increasing  the 
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si/e  of  the  chambers.  Reduction  of  interstage  losses  lead  directly  to  improved  compressor 
and  turbine  overall  efficiencies,  especially  at  higher  operating  pressures. 

Sealing  technology  is  an  area  now  experiencing  increased  importance.  As  the 
operating  pressures  have  increased  in  turbomachinery,  their  tendency  to  leak  has  increased 
since  for  the  same  clearance  level,  a  higher  pressure  ratio  leaks  more  than  a  lower  pressure 
ratio  (until  the  airflow  through  the  seal  becomes  choked).  Some  progress  has  been  made  in 
non-contact  seal  technology,  but  it  has  significantly  lagged  behind  other  technological 
improvements  according  Wrigley  fl|.  In  the  U.S.,  NASA  Lewis  Research  Center  has  also 
pursued  experimental  evaluations  to  improve  sealing  technology,  the  motive  for  this 
research  has  been  reducing  fuel  consumption  of  gas  turbine  engines  for  the  U.S.  Military 
thus  reducing  operating  costs.  According  to  the  NASA  Lewis  report,  contained  in  the 
AGARD  collection  of  1978,  component  eimdencies  still  have  potential  for  improvement, 
and  much  of  that  improvement  could  be  achieved  through  advancements  in  sealing 
technology.  Advanced  seals  could  easily  result  in  a  1/2%  efficiency  gain  for  far  less  cost 
and  risk  than  a  complete  redesign  of  the  compressor  blades  (as  described  above)  1 1  ]. 

In  order  to  improve  sealing  technology,  a  test  bed  is  required  to  better  understand 
the  highly  complicated  flow  through  a  seal,  and  to  adequately  collect  data  in  support  of 
Computer  Fluid  Dynamics  (CFD)  simulations.  To  provide  a  facility  in  support  of 
improving  sealing  technology,  this  two  dimensicmal  frianer  pressurized  air  labyrinth  seal 
test  rig  has  been  designed.  The  major  goals  for  this  test  rig  are: 

1 .  To  coiKluct  experiments  at  Reynolds  and  Mach  numbers  that  are  typical  of  modem 
gas  turtxnes  on  large  scale,  geometrically  similar  2D  seals. 

2.  Support  accurate  bulk  flow  measurements  to  verify  2D  representatiems  of  optimum 
non-contact  seal  designs  through  modification  d*  the  number  of  knives,  roughness, 
orientation,  spacing  and  geometry. 


3.  Make  detailed  flow  Held  measurements  of  sufficient  resolution  by  LDV,  hot  wire, 
and  pressure  transducers  to  compare  with  CFD  code  results. 

This  repext  is  organized  into  five  main  chapters  and  is  followed  by  a  section  of 
riguics  containing  the  test  rig  design  drawings  and  flow  prediction  graphs.  Cha|Xer  I 
provides  an  introduction,  and  is  fdlowed  by  Chapter  1 1 -Background,  which  provides  a 
literature  review  of  labyrinth  seal  research.  Chapter  Ill-Design  Calculations  covoa  design 
constraints,  desired  airflow  inlet  conditions,  test  rig  preliminary  design  proposals,  and 
design  of  each  test  rig  component  including  the  piping  system.  Chapter  III  also  contains 
examine  calculations.  Chapter  IV-Data  Acquisition  System  covers  the  fxoposed  computer 
controlled  data  collection  system  and  its  associated  components  for  which  the  rig  is 
designed  to  emfdoy.  Chapter  V-Test  Rip  Performance  discusses  the  factors  pertaining  to 
labyrinth  seal  effectiveness  and  provides  a  model  formulation  from  which  a  computer 
program  was  written  to  predict  Mach  number,  Reynolds  number,  and  flow  rates  for 
various  seal  configurations.  This  chapter  also  discusses  the  results  of  various  simulations. 
Chapter  VI-CcxiclusicHis  and  Recommendatirxis  describes  the  overall  results  of  this  report 
and  provides  recommendations  for  further  areas  of  study. 
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II.  BACKGROUND 


The  concept  of  limiting  fluid  leakage  using  a  number  of  restrictors  is  not  a  new 
idea.  In  his  review  of  Labyrinth  Seal  literature,  Sneck  [2]  credits  C.A.  Parsons  with 
development  of  the  labyrinth  seal  in  concert  with  Parson's  [3]  development  of  the  steam 
turbine  near  the  end  of  the  1800's.  As  Sneck  describes  it.  Parson's  idea  was  to  create  a 
torturous  flow  path  between  high  and  low  pressure  regions  by  means  of  a  series  of  non¬ 
contacting  restrictors  and  separating  chambers.  The  main  function  of  the  restrictors  was  to 
dissipate  the  kinetic  energy  of  the  fluid  flow  in  the  separating  chambers  and  thus  reduce  the 
leakage  from  high  to  low  pressure  areas.  Parson's  original  method  of  multiple  sealing 
knives  has  changed  very  little  except  for  more  modem  geometry  and  orientations. 

The  first  pioneering  theoretical  paper  on  labyrinth  seal  design  was  introduced  by 
Martin  [4]  in  1908.  Martin  considered  the  labyrinth  to  be  a  series  of  discrete  throttling 
processes  very  similar  to  flow  through  orifices.  Martin  derived  a  formula  for  labyrinth  seal 
leakage  based  upon  a  number  of  simplifying  assumptions  namely:  flow  through  the  seal 
was  isothermal,  all  knives  and  cavities  were  symmetrical,  kinetic  energy  was  completely 
destroyed  in  each  cavity,  and  the  airflow  was  sub  critical  (not  sonic,  i.e.,  flow  not  choked). 
Most  researchers  following  Martin's  published  work  focused  on  modifying  some  or  most 
of  Martin's  original  assumptions.  Gercke  [5]  modified  Marlin’s  equation  with  the  addition 
of  a  kinetic  energy  carry-over  factor.  He  assumed  the  flow  through  the  seal  was  adiabatic  in 
that  the  flow  would  return  to  an  isothermal  condition  via  a  constant  pressure  process  at 
each  throttling.  This  theory  supported  the  establishment  of  his  kinetic  energy  carry-over 
factor.  Gercke  also  modified  Martin's  formula  for  varying  areas  between  the  restrictors. 

The  next  major  modifier  of  Martin's  formula  was  Egli  [6]  in  1935.  Egli  examined 
the  use  of  Marlin's  formula  for  fewer  than  4  restrictors  in  both  incompressible  and  semi- 
compressible  flows.  One  of  Egli's  main  contributions  was  that  he  noted  the  pressure  ratio 
across  each  restrictor  increases  as  the  flow  moves  from  the  first  restrictor  to  tiic  last. 


Therefore  the  last  restrictor  is  always  the  first  (assuming  same  gap  size  through  the  seal)  to 
choke  (reach  cTitical  flow  pressure  ratio).  Egli  also  experimentally  determined  that  the  last 
restrictor  would  not  choke  at  the  same  pressure  ratio  as  a  sharp  edged  orifice  since  its  flow 
coefficient  increases  as  the  pressure  ratio  increases.  Egli  produced  test  results  for  seals 
showing  the  effect  of  gap  size  and  knife  thickness  on  seal  performance. 

Hodkinson  f7|  in  1939  was  the  first  to  approach  the  kinetic  energy  carry-over  factor 
of  Gercke  from  a  fluid  dynamic  standpoint  rather  than  a  thermodynamic  standpoint  He 
considered  the  interaction  of  the  flow  encountering  the  next  sealing  knife  as  a  jet  with 
portions  being  sheared  off  by  the  obstruction  of  the  knife.  Jerie  [8]  in  1948  examined  the 
effect  of  the  knife  thickness  to  clearance  ratio  on  the  seal  performance.  He  determined  that 
when  this  ratio  was  greater  than  2: 1  than  the  knife  behaves  similariy  to  a  rounded  nozzle 
instead  of  a  sharp  edged  orifice.  Jerie  also  determined  that  a  tooth  depth  to  tooth  spacing 
ratio  of  slightly  less  than  1  yields  an  optimum  performing  seal  (minimal  leakage).  There 
were  some  other  papers  between  1948  and  1961  of  minor  importance  from  Kearton  and 
Keh  [9]  and  Zabriske  and  Sterlicht  [10|  coiceming  flow  coefficients  and  friction  factors, 
respectively. 

The  next  major  paper  was  by  Vermes  ( 1 1  ]  in  1961 .  Vermes  crmibines  the  efforts  of 
Martin,  Gercke,  Egli,  and  Zabriske  into  one  major  work.  Vermes  uses  Martin's  formula 
and  derives  new  theoretical  and  semi -theoretical  formulas  fca*  computation  of  the  leakage. 
These  formulas  correlate  within  5%  to  experimental  tests  for  three  types  of  seals  (straight, 
stepped  and  comlxnation).  Vermes  also  calculates  off  (tesign  seal  performance  of  the  seals 
from  both  a  theoretical  and  experimental  standpoint  In  1977  Stocker  (12],  working  for 
NASA  Lewis  Research  Center,  investigated  the  influence  of  rotation  on  seal  leakage  with 
solid-smooth,  abradeable,  and  honeycomb  lands.  He  also  made  experimental 
measurements  of  flow  leakage  for  the  2D  and  3D  cases.  His  conclusions  were  that 
advanced  seal  designs  increase  seal  effectiveness  by  over  25%  when  compared  to 
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conventional  stepped  seals,  and  rotational  effects  produced  mixed  results  with  regard  to 
seal  leakage  depending  upon  the  materials  used  in  constructing  the  seal. 

The  AGARD  conference  of  1978  [1]  specifically  states  that  meaningful  work  in  the 
area  of  sealing  technology  has  not  kept  pace  with  advances  in  the  major  gas  turbine 
component  assemblies,  and  that  resources  need  to  be  dedicated  to  further  woric  in  this  area. 

Chapter  V-Test  Rig  Performance  fully  develops  the  equations  for  determining  the 
flow  through  a  seal  and  describes  the  factors  affecting  seal  leakage,  therefore  these 
equations  will  not  be  developed  here. 
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III.  DESIGN  CALCULATIONS 


A.  CHAPTER  OVERVIEW 

The  test  rig  design  began  with  first  establishing  a  set  of  design  goals.  These  goals 
were  used  to  consider  a  set  of  design  options  and  tradeoffs  and  to  guide  the  establishment 
of  a  set  of  preliminary  designs.  The  preliminary  designs  led  to  selection  of  one  "best 
candidate"  design  for  which  a  more  in-depth  stress  analysis  and  cost  analysis  were 
performed.  After  ensuring  the  design  had  some  promise,  a  complete  stress  analysis  was 
conducted  and  detailed  mechanical  design  produced.  Meanwhile,  a  set  of  performance 
predictions  was  generated  to  verify  the  design  would  achieve  the  original  design  set  points 
with  a  low  level  of  risk  (both  technical  risk  and  safety). 

The  sections  included  in  this  chapter  are  Design  Goals,  Desired  Inlet  Conditions, 
Design  Methodology,  Possible  Design  Approaches,  Air  supply.  Piping  and  Valves,  and 
design  of  the  Diffuser,  Contraction,  Test  and  Exit  sections.  The  chapter  is  completed  by  a 
discussion  of  the  labyrinth  seal,  and  some  design  configurations  used  to  modify  seal 
perfcxmance. 

B.  DESIGN  GOALS 

The  first  choice  to  be  made  was  the  speciHc  design  objectives,  i.e.,  what  type  of 
seals  did  we  desire  to  test  (straight,  stepped,  planer,  annular,  etc.)  and  to  what  extent  did 
we  want  to  test  them.  In  most  seals  the  ratio  of  gap  clearance  to  rotor  diameter  is  so  small 
the  3D  annular  seal  may  be  considered  2D  planer  (Figure  1.1),  especially  if  the  rotor 
rotation  is  not  particularly  high  and  the  fluid  fiiction  coefficient  is  low  (as  is  the  case  with 
air).  With  the  desire  to  match  gap  Reynolds  numbers  and  Mach  numbers  that  occur  in  real 
gas  turbines  (with  realistic  seal  geometry  and  large  scale  to  facilitate  instrumentation)  it  was 
decided  that  an  effective  2D  planer  seal  test  rig  would  suffice.  The  typical  Reynolds  number 
range  in  a  gas  turbine  engine  will  run  to  about  10,000  in  the  compressor  and  to  greater  than 
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100,000  in  the  turbine,  with  axial  Mach  numbers  less  than  0.9  and  pressure  ratios  of  2  to 

3.  A  2D  annular  design  (fixed)  was  considered  due  to  the  lack  of  end  wall  effects  (end  wall 
effects  are  particular  to  2D  test  rigs  due  to  boundary  layer  growth  along  the  edges), 
however  the  ability  to  construct  a  curved  view  port  was  in  doubt  as  was  the  ease  of 
instrumenting  the  rig  (the  rig  would  have  to  be  very  large  to  facilitate  instrumentation),  this 
rig  would  be  cost  |xohibitive.  The  following  majex  design  goals  and  are  listed  below; 

1 .  Safety  was  the  overriding  concern  on  this  project  Thus  a  minimum  Factor  of 
safety  (FOS)  of  3  against  yield  was  the  design  point  for  each  component. 

2.  A  test  rig  of  the  smallest  size  while  preserving  good  aspect  ratio  and  adequate  LDV 
control  volume  was  desired,  this  resulted  in  seals  scaled  upward  by  about  a  10: 1 
margin. 

3.  Continuous  run  operation,  the  airflow  supply  was  limited  by  available  equipment  to 
2000  SCFM  and  150  psig  as  shown  in  Figure  3. 1 .  Blow  down  operation  was 
considered,  however  this  would  impact  other  users  of  the  air  supply  system  by 
depleting  the  reserve  air  supply  and  would  require  recharge  intervals  (however 
users  of  the  blow  down  rig  would  not  significantly  affect  this  rig's  operation  as  it 
would  require  no  recharge  time). 

4.  Cost  and  ease  of  construction  were  very  important,  the  rig  needed  to  be  cost 
effective  and  "consUaction  friendly". 

5.  Once  constructed,  the  test  rig  required  maximum  expandability,  user  friendliness, 
and  low  maintenance.  This  necessitated  the  development  of  a  modular  design. 


8 


C.  DESIRED  INLET  CONDITIONS 


Gas  turbine  engine  labyrinth  seals,  although  performing  the  same  sealing  function 
in  both  the  compressor  and  turbine  sections,  often  encounter  widely  variable  (xmditicms  of 
airflow  conditions  based  upon  design,  operating  speed,  operating  temperature  and 
pressure.  In  fact,  even  seals  within  the  same  section  such  as  two  sets  of  comfniessor  seals, 
may  experience  variable  inlet  air  profiles.  With  this  in  mind,  it  was  desired  for  the  test  rig  to 
possess  some  method  able  to  produce  an  adjustable  air  entrance  profile  at  the  test  secticm 
thus  replicating  real  engine  conditions.  The  effect  of  air  entrance  profiles  on  seal 
performance  is  not  well  understood,  thus  it  is  an  important  aspect  to  consider  in 
constructing  labyrinth  seals. 

The  rig  was  thus  required  to  produce  flow  from  smooth  and  uniform  profiles  to 
turbulent  linear  and  non-uniform  profiles,  with  provisions  for  introducing  shear  or  vortex 
conditions.  A  method  was  also  desired  to  vary  pressure  ratios  from  1.2  to  3.0. 

D.  METHODOLOGY 

To  provide  the  range  of  profiles  described  in  section  C.  above,  a  source  of  airflow 
modification  methods  was  required.  Although  many  books  on  wind  tunnels  were 
reviewed,  only  one  particular  wind  tunnel  book  by  Pope  1 13]  was  consulted  and  its  impact 
on  the  design  of  the  test  rig  appears  in  following  sections.  The  wind  tunnel  information 
was  utilized  to  determine  the  impact  of  screens,  honeycombs,  expansions  and  contractiems 
on  modifying  the  airflow.  Using  the  information  from  the  wind  tunnel  literature  and 
considering  the  constraints  of  part  B,  3  designs  were  considered.  These  designs  were 
evaluated  for  structural  integrity  using  simple  plate  equations  for  stress  and  strain  and 
modified  as  required.  A  final  design  was  chosen  based  upem  strength,  size,  ease  of 
construction,  and  case  of  modification.  Each  of  the  preliminary  designs  is  discussed  in 
.section  E. 
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E.  POSSIBLE  APPROACHES 

Proposal  #1  consisted  of  a  diffuser  designed  for  maximum  pressure  recovery  and 
included  both  horizontal  and  vertical  transitions  of  S-6  degrees  with  internal  splitter  vanes. 
The  diffuser  was  designed  in  2  sections,  a  first  diffuser  stage  followed  by  a  fine  screen, 
then  followed  by  a  second  diffuser  stage  ending  with  multiple  fine  screens  and  a 
honeyoxnb  section.  The  diffuser  was  expected  to  require  35  psig  in  diffuser  assemUy  and 
employed  a  3D  large  area  ratio  (6.7)  final  contraction  section  would  form  the  air  into  a  very 
smooth  profile  (less  than  2%  total  turbulence).  The  problem  with  this  deagn  was  its 
complexity,  projected  construction  time,  size,  cost,  and  possible  problems  with  separated 
or  highly  ncxi-uniform  flow. 

Proposal  #2  was  similar  to  #1  with  the  exception  of  a  modified  diffuser  assembly 
that  would  now  be  one  piece  with  separate  2D  horizontal  and  vertical  transitions,  an  8 
degree  diffusion  angle  with  many  s|4itter  vanes,  and  a  series  of  diffuser  stage  screens  still 
ending  with  1  honeycomb  section.  At  the  diffuser  exit  would  be  a  smaller  contraction  ratio 
(3.3  vice  6.7)  operate  with  35  psig  in  the  diffuser  assembly  (still  designed  for  maximum 
pressure  recovery)  and  maintain  a  3D  final  contraction  section.  The  major  problems  with 
this  design  was  complex  placement  the  numerous  vanes,  still  time  consuming  to 
construct,  and  most  importantly  the  questionable  ability  to  “shear”  the  gauss  distribution  of 
the  air  inlet  into  equal  parts  for  effective  diffusion.  This  design  would  have  even  more 
chance  of  being  affected  by  separated  ot  highly  non-uniform  flow  than  proposal  #1 . 

Proposal  #3  incorporated  a  diffuser  being  of  a  "dump"  type  (not  designed  for 
maximum  pressure  recovery)  with  a  removable  top  plate  for  easy  repositioning  of  screeas 
and  honeycombs  according  to  guidance  received  by  Dr.  Homung  [14].  The  diffuser  would 
contain  only  vertical  transitions  with  no  vanes,  but  a  series  of  torturous  air  paths,  multiple 
diffuser  stage  screens,  and  multifile  honeyccxnb  sections.  Diffuser  operating  pressure 
would  now  be  45  psig  due  to  the  increased  diffuser  losses,  this  was  an  increase  from  35 
psig.  The  final  contraction  section  changed  to  only  a  vertical  contraction  (Area  ratio  of  6). 
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The  problems  with  this  design  are:  it  is  now  very  heavy,  it  must  be  corrosion  protected, 
aiKl  there  is  no  way  to  correctly  predict  the  airflow  prc^ile  and  initial  perforation  plate 
placement.  The  test  section  did  not  change  appreciably  from  the  initial  deagn,  it  is 
constructed  of  steel  and  1"  glass  plate.  Sensors  are  mounted  either  on  seal  or  cm  top  plate 
(replacing  glass),  and  are  discussed  more  fully  in  Chapter  IV. 

F.  AIR  SUPPLY,  PIPING  AND  VALVES 

Air  is  supplied  to  the  test  rig  area  from  a  3  stage  600  Hp  Elliot  centrifugal  air 
compressor  producing  2000  SCFM  at  150  psig  and  180  <>F  maximum  outlet  temperature 
(Figure  3.1).  The  air  passes  through  an  air  drier  assemUy  and  to  an  additional  air 
compressor  capaUe  of  boosting  the  exit  pressure  to  300  psig.  In  most  circumstances  the 
boost  pump  is  bypa.ssed  and  the  air  is  fed  directly  into  an  8000  ft^  storage  tank  group  that 
supplies  the  test  rig  area.  The  air  is  piped  into  the  test  rig  area  through  4  ittch  diameter 
stainless  steel  piping,  the  cormection  designed  to  be  utilized  by  this  test  rig  terminates  with 
a  steel  blank  flange  with  a  standard  6  bolt  pattern  of  1/2"  bolts.  Since  the  supply  manifold 
runs  the  length  of  the  building  and  supplies  other  test  tigs,  the  upstream  building  shutoff 
valve  can  not  be  used  as  the  sde  on/off  valve  for  this  test  rig  since  it  impacts  all  air  users  in 
the  building.  Therefore,  a  dedicated  4"  full  shut  off  valve  rated  to  300  psig  is  required.  The 
shut  off  valve  can  also  be  used  as  a  throttling  valve  when  the  manifold  is  pressurized  to  300 
psig  (when  the  2X  boost  pump  is  in  operation)  to  prevent  pressures  greater  than  150  psig 
in  the  test  rig  piping. 

A  drawing  of  the  pipe  connection  flange,  required  to  replace  the  manifold  blank 
flange,  is  included  as  Figure  3.3.  This  flange  is  designed  for  minimal  flow  loss  and 
reduces  the  4"  manifold  diameter  to  a  2"  diameter.  To  determine  the  flow  loss  coefficient 
for  this  flange  design,  tables  from  While  f  15]  were  consulted.  Thus  for  a  loss  coefficient  of 
0.005,  a  required  radius/overall  diameter  =0.20.  Knowing  the  inlet  diameter  of  the  pipe  is 
1.939",  a  value  of  r  is  determined  to  be  0.388".  The  next  larger  size  (using  1/16" 
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increments)  is  7/16”  that  equals  0.437S”.  Utilizing  the  supplied  airflow  as  a  guide  (2000 
SCFM),  air  mass  flow  was  determined  from  the  following  equation: 


Ptui^'wr  (3-1) 

ft 

(.076^)(2000^)  -  2.55— 
ft  min  sec 

A  pipe  diameter  of  2"  was  chosen  to  reduce  the  cost  of  die  compcments,  and 
facilitate  structural  design  (smaller  structures  are  less  massive  for  equivalent  strength).  In 
the  interest  of  having  a  smooth  pipe  wall  with  minimal  upkeep,  and  simplifiid  pipe 
connections,  the  pipe  chosen  for  this  rig  was  2"  schedule  80  PVC  rated  to  600  psig.  The 
inside  diameter  of  2"  schedule  80  pipe  is  nominally  1.939".  PVC  jnpe  is  very  temperature 
sensitive,  and  loses  strength  with  temperatures  mudi  higher  than  standard  omditions. 
Although  the  com{xessor  oudet  temperature  can  approadi  180  the  massive  size  of  the 
storage  tanks  (8000  ft^  made  of  1.5"  steel  jrfate)  and  long  run  of  piping  exterior  to  the 
building  (apfH^oximately  200  feet)  allows  the  air  supply  manifdd  temperature  to  remain 
below  1 10  op,  especially  with  the  normal  Ntonterey  temperatures  around  65  op.  Schedule 
80  PVC  pipe  must  be  derated  by  33%  strength  at  1 10  op  (to  400  psig),  and  50%  (to  300 
psig)  at  140  op  according  to  ASTM  D-1784. 

The  valves  required  for  this  test  rig  all  must  be  of  the  full  port  design,  i.e.,  the 
inside  air  passage  diameter  must  not  be  less  than  the  pipe  internal  diameter  otherwise  the 
test  rig  air  flow  will  be  reduced  and  the  test  rig  will  fail  to  achieve  designed  operating 
points. 


12 


The  following  valves  (or  equivalent)  should  be  procured  {nrior  to  initial  construction 
in  order  to  account  for  any  changes  required  in  connections  prior  to  construction,  the  valves 
should  be  of  sturdy  construction  and  rated  to  300  psig  (shut  off  and  regulator).  Fisher 
series  667-ET  stainless  steel  valves  would  be  a  good,  but  expensive  choice: 


Name 

Size 

Type 

Notes 

Manual-Shut-off 

valve 

2"  pipe  thread  to  2" 

pipe  thread 

Globe  valve 

fuUptnt 

ability  to  reduce  inlet 

pressure  300  to  150 

Adjustable-Pressure 

RegulatcH- valve 

2"  pipe  thread  to  2" 

pipe  thread 

Regulator 

full  port 

regulate  0  to  45  psig, 

tight  shut  off. 

Pressure-Relief 

valve 

r'pipe  thread 

Relief,  spring 

adjustable 

ability  to  set  relief 

range  35-50  psig 

Outlet-Ctontrol 

valve 

4"  flange  or  4"  pipe 

Throttling  Butterfly 

or  Globe,  full  pon 

max  pressure  30 

psig 

Table  3.1 


A  filter  is  included  in  Figure  3.2,  this  filter  is  required  since  no  upstream  filtration 
capability  exists  within  the  air  supply  system.  The  use  of  the  filter  will  be  required  to 
ensure  long  life  of  all  fixed  sensors,  particularly  hot  wire  probes.  The  filter  chosen  is  a 
Dollinger  coalescing  model  GP- 198-0030-020  rated  for  180  psig  (maximum)  and  840- 
1008  SCFM  at  150  psig  with  a  pressure  drop  of  2  psig.  This  filter  is  advertised  to  have  a 
99.9%  efficiency  at  removing  .3  micron  and  larger  contaminants  (including  water  and  oils). 
A  differential  pressure  gauge  is  included  with  the  filter  assembly  and  indicates  when  the 
filter  media  requires  changeout.  Inlet  and  outlet  connections  are  2"  NPT. 
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G.  DIFFUSER  ASSEMBLY 


To  facilitate  ease  of  construction,  potential  modification,  and  robust  design,  a 
rectangular  dun^  diffuser  is  onployed.  In  cmler  to  adjust  the  air  profile  exiting  the  diffuser, 
a  removable  bolted  top  is  included  (widi  lifting  padeyes).  The  diffuser  can  be  modified  with 
a  number  screens,  honeycombs,  and  perforations  if  desired.  The  factOTS  govoning  the 
application  of  screens,  honeycombs,  perforations,  and  contractitHi  ratios  are  discussed  in 
section  H.  Diffuser  section  drawings  ate  contained  in  Figures  3.4  to  3.10. 

The  Diffuser  assembly  is  comprised  of  a  1"  steel  bottom  plate  with  0.75"  welded 
side  plates  and  a  0.75"  welded  entrance  plate.  These  plates  have  welded  flange  connections 
to  secure  the  removable  1"  top  plate.  Internally  the  diffuser  holds  welded  sections  of  0.25" 
steel  plate  of  various  lengths  to  act  as  spacers  for  the  perforation  plates,  screens,  and 
honeycombs.  Allowance  is  made  for  two  side  plates  and  a  bottom  plate  welded  together  to 
form  each  spacer,  the  exit  section  flange  is  stepped  down  to  secure  the  last  set  of  spacer 
plates.  The  exit  of  the  diffuser  section  has  a  bolted  connection  flange  to  the  contraction 
section.  All  bolts  used  in  the  diffuser  and  contraction  sectitm  are  1/2"  diameter  SAE  grade  1 
or  above,  coarse  threaded.  The  following  factcx^  of  safety  (FOS)  exist: 

3.3  diffuser  front  plate  (FOS  9.9)  with  3X  stress  concentration 

factor. 

4.42  diffuser  top  plate 
9.9  diffuser  side  plates 
5.75  attachment  bolts 

In  completing  stress  calculations  for  the  top  plate,  an  expected  maximum  operating 
pressure  in  the  diffuser  was  chosen  as  45  psig.  The  diffuser  top  plate  area  subject  to  this 
internal  air  pressure  was  an  area  18"  wide  and  33"  in  length.  Using  Roark's  Handbook  for 
Stress  and  Strain  [16]  in  the  case  of  a  rectangular  plate  with  4  fixed  edges  (case  8)  the 
following  formulas  apply; 
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R)r  maximum  stress  at  the  cotter  of  the  long  edge: 


0.4974yfc^ 


(3-2) 


For  maximum  stress  at  the  center  of  the  plate: 

0.2472<?f»* 

- — ;3— 

For  maximum  deflection  at  the  center  the  [date: 

().02»qb* 


Variables:  q  -  45psig 

E  -  30E6psi 
t  ■  thickness 


(3-3) 


(3^) 


From  these  formulas,  the  maximum  stress  at  45  psig  on  the  1  ”  top  plate  is  7246  psi 
at  the  edge,  3601  psi  at  the  center,  with  a  maximum  center  deflection  of  0.01”.  The  same 
formulas  are  utilized  for  the  surrounding  [dates  to  generate  stress  predictions.  To  calculate 
b(dt  strength,  Shigley's  Mechanical  Engineering  Design  book  ( 17]  was  utilized  for  bolt 
specifleations.  Shigley's  book  contains  specific  formulas  to  optimize  the  use  of  bolts  and 
fasteners  based  upcm  their  stiffness,  however  his  method  was  very  time  consuming  and 
cumbersome.  Through  previous  experience  with  these  formulas  it  was  determined  that  by 
choosing  a  significant  factor  of  safety  (generally  4  or  more)  against  the  proof  yield  strength 
of  the  bolts,  a  very  conservative  and  yet  rapid  determination  cd*  the  number  of  fasteners 


could  be  made.  With  this  in  mind,  an  initial  FOS  of  5  was  chosen,  and  based  upon 
geometric  considerations  the  resulting  number  of  1/2"  bolts  (using  the  lowest  strengdi  S  AE 
grade  1)  was  determined  to  be  34,  this  in  turn  restilted  in  the  5.75  FOS. 

Force  on  top  plate=(45  psi)(18")(33")=26,730  Ibf 
1/2"  bolt  stress  areasO.142  in^ 

Grade  1  yield  strengdi=32  kpsi 
(0.142  in2)(32kpsi)=4544  Ib^lt 
so  for  FOS=5,  each  bolt  can  only  hold  909  Ibf 
therefore  #  bolts  =26730/909=29.4  bolts 

Screens  and  honeycombs  are  the  primary  flow  conditioners  used  to  noodify  airflow 
in  the  diffuser,  and  they  may  be  arranged  at  various  positions  within  the  diffuser  assembly 
utilizing  0.25"  spacer  plates  as  discussed  previously.  Although  the  perforated  sections  will 
aid  in  breaking  up  the  Gaussian  air  flow  {nofile  exiting  the  inlet  pi{)e  to  the  diffuser,  no 
reliable  data  exists  concerning  the  best  configuration  fen  this  task.  The  {neferred  location  of 
the  1st  {)erforation  plate  is  even  with  the  inlet  tube  cap  at  approximately  5.0"  from  datum. 
Datum  is  the  zero  {wint  defined  as  the  inlet  wall  of  the  diffuser.  The  perforated  inlet  pipe 
attempts  to  double  the  area  of  the  mean  inlet  flow  and  strain  in  to  the  sides  of  the  diffuser. 
The  1st  {)erf(nation  plate  has  twice  the  airflow  area  of  the  {noceeding  {lerforated  section  and 
tries  to  maintain  the  flow  towards  the  edges  of  the  diffuser.  It  is  located  at  the  9.25"  datum 
location.  The  second  i)erforation  plate  again  doubles  the  effective  flow  area  and  tries  to 
form  a  somewhat  uniform  profile  flow  into  the  honeycomb  section.  The  second  i)erforated 
plate  is  located  at  9.5"  from  datum. 
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Screens  have  seen  very  widespread  use  both  as  turbulence  reducers  (fine  wire  and 
spacing  far  upstream  of  the  test  section)  and  as  turbulence  enhancers  (coarse  wire  cloth  just 
preceding  the  test  secticm).  Screens  typically  behave  in  the  relation  as  set  forth  by  Dryden  in 
1929  [18]  ; 


U'^  _  1 

U,  (1 + 


(3-5) 


where  (/],and  U\^  are  the  root  mean  square  of  the  velocity  fluctuations  downstream  and 
upstream  of  the  screen,  respectively  in  x,  y,  and  z  directions.  The  pressure  drop  factor  k, 
and  the  number  of  screens  n  are  also  variables.  As  an  example  problem,  if  2  screens  are 
used  and  each  screen  is  comprised  of  18  mesh  wire  of  0.011"  diameter  wire  with  a 
k=0.842  at  a  velocity  of  20  feel/second,  the  turbulence  will  be  reduced  by  46.7%. 

A  concern  when  utilizing  screens  is  the  drag  on  the  mesh  due  to  the  airflow.  In  the 
above  case  the  drag/ft^  =kq,  where  q  is  the  mass  flow/ft2,  the  drag  is  0.41b/ft2  and  for  a 
screen  of  9"X18"  the  total  drag  would  be  only  1.13  Ibf  on  the  screen.  The  recommended 
screens  for  the  test  rig  should  be  constructed  of  20  mesh  0.017"  wire  diameter  (with 
K=1.8)  screen  located  at  the  20",  24",  and  28"  points  from  diffuser  datum.  This  spacing 
will  allow  adequate  mixing  and  evening  out  of  the  flow  between  the  screens.  These  screens 
should  yield  a  tiffbulence  reduction  of  78.7%  using  Equation  3-5. 

Honeycombs  are  primarily  used  to  strain  the  flow  and  make  it  close  to  ID  as 
possible  by  breaking  up  large  scale  structures.  TTie  honeycomb  prevents  circulation  and 
lateral  velocity  variations  from  propagating  as  the  flow  progresses.  The  only  real  rules  of 
thumb  for  the  honeycomb,  is  that  the  length  of  the  honeycomb  should  be  5  to  10  times  its 
cell  diameter,  and  the  octagon  configuration  typically  yields  the  lowest  pressure  drops 
compared  to  circular  or  square  patterns.  The  recommended  honeycomb  is  0.25"  diameter 
hexagons  of  aluminum  construction  of  2.5"  in  length.  The  honeycomb  leading  edge  should 
be  located  at  the  13.5"  datum  location  and  secured  by  press  fitting  into  place  and  placing 
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0.25”  spacer  plates  upstream  and  downstream  of  this  location.  Wall  interfaces  should  be 
sealed  by  the  use  of  a  flexible  RTV  or  Permatex  to  prevent  errant  flow  around  the 
honeycomb. 


H.  CONTRACTION  SECTION 

The  contraction  ratio  of  any  test  rig,  or  wind  tunnel,  is  an  important  item  to  consider 
for  two  reasons:  1)  the  desired  contraction  ratio  will  determine  the  size  of  the  diffuser  and 
test  sections,  and  2)  the  contraction  ratio  has  a  major  effect  on  the  longitudinal  velocity 
variations  within  the  flow  path.  As  discussed  by  Pope  [13],  a  properly  chosen  contraction 
ratio  has  a  major  contribution  to  improving  the  performance  of  a  wind  tunneiytest  rig.  To 
determine  the  effect  of  contracdon  ratio  some  background  is  required  in  accordance  with 
Pope's  discussion. 

Vu=upstream  velocity  variation 
vd=velocity  variation  after  the  contraction 
Vu=upstream  mean  velocity 
Vd=mcan  velocity  after  the  contraction 
n=contraction  ratio 

Writing  Bernoulli's  equation  assuming  pressure  and  gravity  loss  is  negligible 

iP(V.  +  v.)"=|p(V',  +  v,)^  (3-6) 


v>2v.V.=vj  +  2v,V, 


(3-7) 


Dividing  by  Vl  /  Vl  and  neglecting  the  (v/V)2  term  as  being  very  small 


2k. 


(3-8) 


18 


(3-9) 


V*/ Vj  -1//I* 


Therefore,  the  variation  ratio  of  upstream  velocities  are  variable  through  the  inverse 
of  the  square  of  the  contraction  ratio.  For  a  contraction  ratio  of  6,  the  mean  velocity 
variaticMis  would  be  reduced  to  1/36  o(  their  initial  value,  or  yield  a  97.2%  reducticxi.  The 
contraction  section  profile  was  chosen  from  Pope's  book  on  Wind  Tunnel  Oesign  fl3]. 
Pope  also  recommends  the  axitraction  cone  be  led  by  a  settling  chamber  equal  to  0.5  of  the 
equivalent  air  chamber  diameter,  and  followed  by  settling  chamber  of  the  same  size  just 
pritx  to  the  test  section.  A  high  quality  inside  oxitiaction  profile  graph  is  included  with 
coordinates  as  Figure  3.12,  this  figure  should  be  used  to  shape  the  ccmtiaction  curvature 
with  a  9.0”  distance  between  top  and  bottom  plates  at  the  inlet  and  a  distance  d*  l.S” 
between  the  same  plates  at  the  outlet 

I.  SETTLING/MODIFICATION  SECTION 

This  portion  d*  the  test  rig  allows  the  air  to  "rest”  and  restabilize  after  contracting  in 
accordance  with  Pope's  discussion  in  fI3|.  For  this  aze  of  airflow  area  (1.5”X18”)  an 
effective  length  of  2.94”  is  recommended.  The  actual  settling  camber  available  length  is  5” 
with  an  additional  1”  provided  prior  to  the  test  section  view  area.  This  3”  of  available  space 
can  be  used  to  generate  specific  profiles  as  desired  by  the  user.  Some  specific 
recommended  flow  patterns  zuid  how  they  can  be  generated  are: 

1 .  Shear  [xotile  -  achieved  by  placing  a  fine  screen  covering  only  a  portion  of  the 

airflow  area  and  leaving  the  remaining  area  free  from  obstructions. 
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2.  Vortex  pfoTile  -  using  a  series  of  partial  qjan  win^ets  to  generate  any  number  of 
wingtip  vortices  (generating  axial  vortex  cores). 

3.  Turbulent  profile  -  achieved  through  the  insertion  d*  a  series  of  trip  wires  to  create 
vortex  shedding  and  turbulence,  or  turbulence  generated  by  the  injection  of  air 
normal  to  the  flow  path. 

Example  drawings  of  each  profile  is  contained  in  Appendix  C:  Row  Modifications. 
The  above  recommendations  are  only  a  few  of  the  many  different  profiles  which  could  be 
achieved  not  only  through  modification  of  flow  in  the  settling  section,  but  by  modification 
of  the  flow  pattern  upstream  in  the  diffuser  and  contraction  sections. 

J.  TEST  SECTION 

The  test  section  is  comprised  of  four  major  compements,  two  symmetric  side  plates, 
a  bottom  {^ate,  and  a  top  plate  each  bolted  to  one  another  and  belted  to  the  settling  and  exit 
sections.  The  side  plates  can  be  either  machined  or  welded  out  of  O.S*  steel  plate,  along 
with  the  top  and  bottom  plates.  The  recommended  construction  technique  for  the  side  plates 
however,  is  machining  out  of  1.5”  mild  steel,  with  the  botttxn  and  top  i^ates  being  of 
welded  construction.  The  test  secticxi  allows  for  designed  gaps  of  0. 15"  with  seals  6"  long 
and  18"  wide.  The  test  section  is  designed  to  be  positive  locking  with  no  adjustment  of  the 
seal  once  installed,  the  seal  will  be  bdted  into  place  from  underneath  the  bottom  section. 
Three  glass  top  plates  allows  the  use  of  a  traversing  LDV  system  to  measure  flow  fields 
from  full  front  to  back  of  the  seal  an  virtually  all  area  covering  the  width  of  the  seal.  Each 
glass  plate  is  1  "  thick  with  top  dimensiems  of  6”X4”,  and  bottom  dimensions  erf*  5.5"X4", 
the  fremt  and  rear  surfaces  of  the  plate  are  beveled  to  support  placement  of  the  {rfates  in  the 
test  section  top  plate.  Each  glass  |:4ate  should  be  fx-otected  frexn  direct  glass  to  metal  contact 
by  the  insertion  of  thin  sealing  gaskets  or  putty.  Removal  of  the  top  plate  assembly  fa- 
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cleaning  or  access  to  the  interior  erf*  the  test  section  will  be  accomplished  by  removing  the 
twenty-two  0.375”  bolts.  While  the  top  plate  is  removed  the  glass  plates  should  be  adjusted 
to  give  a  smooth  surface  from  the  front  to  back  of  the  top  plate.  Some  f(xm  of  filler  (RTV) 
should  be  used  to  fill  the  small  gaps  between  the  glass  and  steel  plates.  The  experimental 
seal  is  axistructed  of  steel  to  ensure  robustness  of  the  edges  and  ccxners.  It  can  be  made  in 
1  piece,  or  multiple  pieces  depending  on  the  desires  of  the  user.  The  seal  is  secured  by 
drilling  six  0.025”  holes  in  the  bottom  plate  and  passing  bdts  through  the  bottom  plate  into 
the  seal  that  itself  has  been  drilled  and  tapped.  Additiorud  hedes  for  the  passage  of 
instrumentation  should  be  avoided  without  additional  stress  factor  calculations. 
Instrumentation  cables  should  be  passed  through  to  the  exit  section  where  they  should  exit 
the  rig  at  one  of  the  1  ”  back  pressure  holes  and  be  sealed.  Details  of  the  test  sectiem  design 
are  contained  in  Figures  3. 14  to  3. 16.  The  fcdlowing  factors  of  safety  (FOS)  exist: 


9.9  Flanges  inlet  and  exit 

3.8  Bottom  plate 

6.2  Each  plate  glass 

3.5  Top  plate  steel 

12.5  Side  plates 

6.75  3/8"  attachment  bolts 


K.  EXIT  SECTION 

The  exit  section  is  designed  to  minimize  back  pressure  when  exit  to  atmospheric 
conditions  is  desired,  yet  allow  back  pressure  to  30  psig  as  requited  for  experimentation. 
Details  of  this  design  ate  included  as  Figures  3.17  and  3.18.  For  operation  with  minimal 
back  pressure,  4  additional  1"  diameter  tapped  holes  are  provided  to  vent  internal  pressure 
if  required  and  to  allow  passage  of  the  seal  instrumentation  cables.  These  holes  will  be 
plugged  for  normal  operations. 
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IV.  DATA  ACQUISITION  SYSTEM 


A.  ACQUISITION  SYSTEM  OVERVIEW 

To  effectively  collect,  process,  store  and  display  large  amounts  of  data  from 
multiple  sensors  an  automated  data  acquisition  system  was  required.  NFS  has  recently 
obtained  a  number  of  IBM-PC  compatible  DATASTOR  486-DX  computers  with  14" 
color  monitors  containing  16  megabytes  of  random  access  memory  and  535  megabyte 
hard  drives.  These  computers  have  full  tower  cases  and  six  free  expansion  slots.  The 
school  also  has  numerous  National  Instruments  AT  MIO-16F-5  Multi-function  I/O 
boards,  IEEE  488  General  Purpose  Interface  Bus  (GPIB)  boards,  and  SC2070  GPIB 
termination  breadboards.  The  computers  are  configured  with  Microsoft  Windows  V3.1 
and  run  Lab  view  for  Windows  software.  Labview  is  a  full  featured  graphical  programming 
system  for  data  acquisition  and  control,  data  analysis,  and  data  presentation.  Labview  is 
easily  programmed  to  accept  inputs  from  plug-in  boards  and  supports  numerous  data 
collection  options  beyond  the  scope  of  this  discussion.  This  basic  system  will  be  interfaced 
with  Scanivalve  Corporation  self-calibrating  fixed  temperature  and  pressure  transducers  to 
provided  basic  data  collection  as  described  in  the  following  sections.  This  will  provide  an 
excellent  data  collection/analysis  system.  A  schematic  diagram  of  the  system  is  included  as 
Figure  4.1. 

B.  AIRFLOW  MEASUREMENT 

To  accurately  measure  the  flow  through  the  test  rig,  a  provision  for  a  replaceable 
orifice  plate  has  been  made.  Orifice  units  typically  contain  both  upstream  and  downstream 
pressure  taps  and  utilize  a  calibrated  sonic  orifice  plates.  The  plates  are  replaceable  to 
maintain  accurate  flow  measurement  through  the  pressure  range  of  the  test  rig.  Orifice 
plates  should  generally  be  chosen  to  provide  less  than  1%  deviation  in  the  intended  realm 
of  rig  operation  pressures.  The  upstream  and  downstream  pipe  connections  for  the  orifice 


22 


assembly  are  2"  NPT  of  1.939"  inner  diameter  on  this  test  rig.  The  orifice  assembly  is 
designed  to  be  mounted  just  downstream  of  the  pressure  regulator  valve,  and  just  prior  to 
the  entrance  of  the  diffuser  assembly.  In  this  position  the  orifice  assembly  should 
experience  pressure  less  than  50  psig;  however,  a  potential  drawback  to  this  location  is  the 
high  pressure  drop  induced  in  the  orifice  which  may  prevent  the  test  section  from  achieving 
the  desired  pressure  and  flow  rate.  An  alternative  position  is  to  mount  the  orifice  assembly 
upstream  of  the  pressure  regulator  (150  psig  side),  in  this  position  the  orifice  may  be 
slighdy  less  accurate  due  to  the  larger  operating  pressures,  but  will  guarantee  sufficient  air 
inlet  to  the  pressure  regulator  and  test  section.  The  diffuser  assembly  and  pressure  regulator 
threaded  connections  can  therefore  be  easily  switched  to  allow  placement  of  the  orifice 
assembly  at  the  location  desired  while  preserving  sufficient  upstream  and  downstream  pipe 
lengths  to  yield  good  orifice  results.  A  straight  pipe  length  of  at  least  39"  should  precede 
the  orifice  assembly  for  a  fully  developed  turbulent  air  profile  to  occur,  and  a  trailing 
straight  pipe  length  of  at  least  12"  downstream  of  the  orifice.  Channels  lA&B  through 
4A&B  for  both  temperature  and  pressure  are  dedicated  to  collecting  this  data. 

C.  FIXED  SENSORS 

To  determine  the  condition  of  the  air  at  any  point  in  the  test  rig,  a  series  of  pressure 
taps  and  thermocouples  and  a  traversing  hot  wire  are  employed.  Two  computer  expansion 
slots  are  filled  with  SDIU  MK5  digital  cards,  with  each  card  connected  to  its  associated 
scanivalve.  Although  the  scanivalx  ss  are  driven  by  the  same  drive  system,  one  48  channel 
valve  set-up  is  configured  for  pressure  measurements  (channel  A),  while  the  other  is 
configured  for  temperature  measurements  (channel  B).  Channels  5A  through  8B  are 
dedicated  to  the  diffuser  section,  channels  47A&B  and  48A&B  to  the  exit  section,  and  the 
remaining  38A&B  channels  to  the  test  section.  A  major  design  deficiency  of  the  test  rig  is 
its  current  inability  to  support  rapid  repositioning  of  the  sensors  while  the  unit  is  in 
operation  (however,  most  test  rigs  do  not  have  this  ability).  In  fact,  no  repositioning  is 


23 


currently  possible  without  partial  disassembly  of  the  test  rig.  Sensors  can  be  mounted  at  the 
discretion  of  the  user  depending  upon  the  goal  of  the  experiment,  but  a  recommended 
layout  is  shown  in  Figure  4.1.  In  the  recommended  layout,  the  left  side  of  the  test  rig 
contains  more  in-line  sensors  than  the  right  side,  this  will  allow  the  user  to  determine  if  any 
ill  effects  occur  due  to  this  placement.  Every  data  location  is  also  supported  by  at  least  one 
and  sometimes  two  redundant  sensors.  These  redundant  sensors  should  also  ensure 
adequate  safety  monitoring  capability.  The  scanivalve  system  is  self  calibrating  and  has  an 
advertised  accuracy  of  less  than  0.1%  of  full  scale. 

A  DANTEC  traversing  hot  wire  system  is  also  recommended  to  aid  in  determining 
flow  fields,  and  as  a  supplemental  mass  flow  indicator.  This  system  requires  two  computer 
slots,  one  for  a  GPIB  timing  board,  the  other  for  an  AT  MIO-16F  Analog  to  Digital  board. 
A  calibration  unit  is  provided  to  facilitate  easy  verification  of  data.  This  system  is  software 
controlled  and  is  fully  position  able  through  remote  means  in  three  dimensions.  The  results 
of  using  the  hot  wire  should  be  compared  with  the  LDV  system  described  in  part  C. 
Accuracy  is  advertised  as  less  than  0.1%  of  full  scale. 

D.  NON-INTRUSIVE  LDV  MEASUREMENT 

Many  methods  exist  for  determining  the  state  of  a  flow  Held,  three  of  them  are  hot 
wires,  pitot  probes/pressure  taps,  and  Laser  Doppler  Velocimeters.  Of  the  three  only  one,  a 
traversing  LDV,  can  generate  high  accuracy  data  exterior  to  the  test  rig.  The  NFS's  recent 
acquisition  of  a  dedicated  LDV  system  complete  with  seeding  atomizer  promises  to  yield 
signiHcant  results  when  used  in  concert  with  the  NFS  high  speed  and  low  speed  test  rigs. 
This  LDV  system  is  a  DANTEC  system  with  an  Enhanced  Burst  Spectral  Analyzer,  3 
watt  output  wtaer  cooled  laser,  and  a  stand  alone  data  acquisition  system.  Software  is 
provided  by  the  manufacturer  to  collect,  display,  analyze,  and  modify  the  LDV  data. 

The  LDV  works  on  the  principle  of  Doppler  shifting  to  determine  the  ID  velocity 
of  the  flow  field.  The  laser  beam  is  produced,  then  separated  into  multiple  colored  beams 
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(typically  split  by  a  prism  with  the  two  strongest  beams  of  blue  and  green  collected  at  the 
output).  Each  of  these  beams  are  then  split  again,  with  1/2  of  the  beam  frequency  shifted 
slightly  out  of  phase.  Each  beam  is  then  focused  into  the  control  volume  and  scattered  by 
the  motion  of  the  flow  field,  sometimes  seeding  material  is  required  to  aid  in  scattering  of 
the  beams.  The  beams  are  then  collected  in  either  a  forward  scatter  or  back  scatter  mode 
and  analyzed  to  determine  the  resulting  frequency  shifts.  Forward  scatter  is  typically  the 
better  mode  when  possible,  but  the  test  rig  will  be  limited  to  back  scatter  operation  only. 
These  frequency  shifts  are  then  converted  into  velocities  and  collected  to  obtain  a  velocity 
flow  field. 
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V.  PREDICTED  TEST  RIG  PERFORMANCE 


A.  EQUATION  DEVELOPMENT 

To  accurately  predict  the  leakage  rates  for  a  given  labyrinth  seal  geometry,  and  to 
predict  Mach  numbers  and  Reyncdds  numbers,  a  method  was  required  to  ^ecdvely  predict 
flow  rates  over  large  variations  in  seal  gecxn^ry  and  inlet  pressures.  In  order  to  accomplish 
this,  a  program  was  written  in  MATLAB  (appendix  A)  to  satisfy  these  requirements. 
Section  B  discusses  reduction  of  the  data  into  Madi  and  Reynolds  numbers,  and  the 
speciHc  program  predictions.  In  order  to  develop  the  flow  prediction  program  certain 
assumptions  and  simi^ifications  were  required: 

1.  All  inlet  and  outlet  conditions  of  the  air  are  constant,  known  quantities. 

2.  The  working  fluid  (air)  is  a  perfect  gas. 

3.  Row  through  the  seal  is  adiabatic. 

4.  The  air  flow  behaves  in  a  semi-compressible  fashion  at  speeds  below  Mach  0.9. 

(i.e.  the  calculated  air  density  is  tte  average  of  the  inlet  and  exit  air  densities  across 
the  seal) 

5.  All  jet  kinetic  energy  is  dissipated  at  each  intermediate  seal  cavity. 

6.  The  kinetic  energy  carryover  factor  (P)  and  vena-contracta  modification  (Cc)  are 
constant  for  each  knife  (p=CcP-0.65). 

And  although  not  necessary,  the  program  is  written  with  two  fixed  parameters: 
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I .  The  spacing  between  the  sealing  knives  are  constant 


2.  The  gap  clearance  of  each  knife  is  constant 


The  simultaneous  flow  equations  for  both  ideal  incompressible  and  ideal  soni- 
compressible  flow  will  be  developed  for  one  sealing  knife.  Utilizing  the  mass  flow 
equation  for  ID  steady  flow: 

m^pAV  (5-1) 

and  the  equation  for  the  density  of  air 


P- 


P 

RT 


(5-2) 


and  the  Bernoulli  equation  neglecting  changes  in  height  with  work  terms  set  to  7jen\  and 
treating  the  inlet  velocity  is  zero  due  to  a  sufTiciently  large  reservoir  (denoting  coiKlitions 
upstream  of  the  first  knife  by  the  subscript  ”u"  and  conditions  at  the  downstream  exit  by 
"d"): 

(5-3) 


By  rewriting  terms  Equation  (5-3)  becomes: 


(54) 


Combining  equations  (54)  and  (5-2)  into  Eqiuuion  (5-1)  an  equation  for  the  mass  flow 
through  the  seal  is  obtained: 


To  calculate  the  semi-ixunpressibie  How  across  one  sealing  knife  some 
modifications  to  Equations  (5-1)  through  (5-5)  will  result  Utilizing  the  mass  flow  equation 
(where  the  represents  an  average  of  the  upstream  and  downstream  conditions); 


m*^p*AV* 


(5^) 


and  the  equation  for  the  density  of  air  (assuming  the  air  density  is  now  based  on  the 
average  pressure  across  the  seal): 


(5-7) 


By  rewriting  terms  ot  the  Bernoulli  equation  using  the  values  while  still  denoting 
conditions  upstream  of  the  first  knife  by  the  subscript  "u"  and  conditions  at  the  downstream 
exit  by  "d"  Equation  (5-3)  becomes: 


Combining  Equations  (5-6)  and  (5-7)  into  Equation  (5-8)  an  equation  for  the  semi- 
cxHupressible  mass  flow  through  the  seal  is  obtained; 


m*  -  A 


kP^-Pd)iPu+pj)  A 

I  RT  "  4Wf 


{Pl-Py 


The  ideal  equation  (5-9)  must  now  be  modified  by  a  flow  coefficient  p,  which  is  the 
actual  flow  through  the  seal  divided  by  the  ideal  mass  flow. 


28 


An  expresaon  for  (a  is  the  product  the  contraction  coefficient  Cc  (which  is  the 
ratio  d*  the  minimum  jet  Dow  area  occurring  at  the  vena-contracta  to  the  minimum  area  of 
the  orifice)  and  the  kinetic  energy  carry-over  factor,  p  in  1 17]: 

M^CcP  (5-10) 

The  theoretical  value  of  Cc  is  0.61 1,  but  has  been  experimentally  determined  to  be  0.65  fcH* 
a  sharp  edged  orifice  by  Lamb  ]  19].  An  empirical  relatirai  for  p  based  upon  seal  geometry 
has  been  determined  by  Vermes  [  1 1  ]  as: 


1 

(1-a)" 


(5-11) 


8.52 

- 

-  +  7.23 
d 


(5-12) 


where  I  is  the  distance  between  seals,  and  6  is  the  gap  between  the  top  of  the  knife  and  the 
stator  with  both  surfaces  static. 

In  combining  Equation  (5-9)  with  (5-10)  a  modified  flow  rate  equation  can  now  be 
represented  with  q  as  the  mass  flow  per  unit  width  ot  the  seal. 

«  -  -  py  (5-13) 

In  reality,  the  p  modification  should  be  0.65  fcr  the  first  knife  since  no  kinetic  carry-over 
exists,  with  the  product  of  0.65  and  p  fen-  the  remaining  knives.  To  simplify  equatiem 
development  and  MATLAB  programming,  a  constant  value  of  0.65  is  used  for  all  values 
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of  |i  althcxigh  the  program  does  calculate  attd  display  a  "correct*  value  for  p  based  upon 
Equations  (5-11)  and  (5-12). 

As  displayed  in  Figure  5. 1  (a  four  knife  seal).  Equation  (5-13)  can  be  modified  for 
each  seal  using  the  upstream  and  downstream  pressures.  When  setting  each  q  equation 
equal  to  one-another  and  canceling  the  common  terms  (R,T,p,  and  &  )  the  following 
relationship  is  achieved. 


(Pu2.  p,2)  =  (p,2.  P22)  =  (P22.  P32)  =  (P32.  pd2)  (5.14) 

Equation  (5- 14)  can  now  be  rewritten  into  equations  for  PI ,  P2,  and  P3  in  terms  of  only 
Pu  and  Pd: 


(5-15) 


(5-16) 


(5-17) 


Similar  relationships  can  now  be  determined  for  N  sealing  knives.  Since  the  pressure  ratios 
across  each  knife  can  now  be  determined,  the  use  of  the  compresvsible  flow  equation  for 
isentropic  How  fmm  White  (15]  reveals  the  Mach  number  across  each  seal: 


(5-18) 


Using  the  relationship  Ma=V/a ,  where  a=(YRT)-^  (Y=1  -4  for  air),  the  velocity  of  the  air 

How  through  each  knife  can  be  determined. 
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Reynolds  numbers  can  be  determined  by  using  the  gap  6.  the  velocity  V,  and 
kinematic  viscosity  of  air  v  (with  v  averaged  for  the  pressure  variations  across  the  seal) 
according  to  the  formula 


fi  V 


(5-19) 


B.  PERFORMANCE  PREDICTIONS 

The  test  rig  is  designed  to  operate  with  inlet  pressures  as  high  as  45  psia  and  outlet 
pressures  between  15  psia  and  35  psia  with  seals  of  2  to  5  knives  and  Mach  numbers 
below  0.9.  The  program  in  Appendix  A  was  validated  by  several  methods  for  absolute 
pressure  ratios  to  5  (with  atmospheric  pressure  exit),  and  for  seals  with  between  2  and  14 
knives.  First ,  the  generated  mass  flow  rates  at  the  inlet,  outlet  and  across  each  knife  were 
checked  to  ensure  consistency  (they  agree  to  8  decimal  places).  Second,  intermediate 
pressures  were  reviewed  and  Mach  numbers  checked  to  ensure  increasing  Mach  numbers 
towards  the  downstream  seals.  Third,  the  program  was  run  and  the  output  compared  with 
the  curves  generated  by  the  well  known  Martin's  seal  formula  14]  (modified  by  the  flow 
coefficient  p-0.65  ),  in  which  the  generated  curves  were  virtually  exact  Fourth,  the 
program  was  run  for  a  g^  clearance  of  0.040  inches  and  5  knives  and  the  flow  rate  verses 
pressure  ratio  [dotted.  This  graph  was  compared  to  a  graph  previously  generated  by 
Vermes'  (11]  in  his  work  comparing  a  modified  version  d*  Martins  formula  to 
experimentally  derived  data  Again,  the  MATLAB  program  and  Martin's  fcxmula  were 
virtually  exact  and  correlated  extremely  well  to  Verme's  experimental  data  (Figure  5.2). 
The  MATLAB  program  was  again  run  for  a  seal  with  14  knives  and  a  pressure  ratio  cd*  5  to 
determine  the  intermediate  chamber  pressures  within  the  seal.  These  were  plotted  on  the 
graph  of  pressure  ratio  verses  scaling  knife  number  and  compared  with  a  prior  graph  of 
experimental  data  (again  a  product  of  Vermes'  work)  fll).  The  MATLAB  program 


31 


predicted  slightly  lower  chamber  pressures  than  the  experimental  results  and  hence  a  lower 
rate  d*  seal  leakage  (which  is  to  be  expected  due  to  the  assumptions  written  into  the 
program),  but  the  predicted  curve  differed  from  the  measured  curve  by  only  3%  to  5% 
(Figure  5.3). 

Confident  that  the  MATLAB  program  is  reasonably  accurate  to  within  5%, 
Reynolds  number  and  Mach  number  verses  pressure  ratio  curves  were  generated  for  2,  3, 
4,  and  5  knife  seals  each  with  gaps  of  1/32",  2/32",  3/32",  4/32",  5/32".  A  series  of  four 
graphs  were  produced  with  exit  pressures  of  5  psia,  15  psia,  25  psia,  and  35  psia  (Figures 
5.4,  5.5,  5.6,  and  5.7  respectively).  Figure  5.8  is  actually  Figure  5.4  disf^aying  only  the 
maximum  and  minimum  Reyndds  number  and  Mach  number  ranges.  Figure  5.9  displays 
low  speed  predicticxis  (pressure  ratios  less  than  1.2)  for  an  exit  pressure  of  15  psia. 

Figure  5.10  shows  the  effect  of  varying  back  pressure  on  the  de-coupling  of 
Reyndds  numbers  and  Mach  numbers.  Thus  for  a  given  Mach  number,  a  variation  in 
Reyndds  number  can  be  obtained  by  adjustment  of  the  test  rig  outlet  contid  valve  to 
produce  a  back  pressure  (above  atmospheric  pressure)  on  the  exit  ade  of  the  seal. 
Conversely,  if  an  attachment  were  connected  at  the  exit  to  the  test  section  capable  of 
drawing  a  vacuum,  the  Reyndds  number  range  could  be  further  increased. 

In  summary,  based  on  Figure  5.8,  the  test  section  should  operate  with  satisfaction 
below  Mach  0.9  with  Reynolds  numbers  up  to  1(X),000.  Figures  5.11  through  5.14 
display  the  intermediate  chamber  pressures  and  seal  leakage  rates  for  23,4,  and  5  sealing 
knives  as  a  function  d  pressure  ratio  (to  3)  and  atmospheric  exit  pressure.  Those  figures 
graphically  show  that  Mach  number  is  always  highest  across  the  last  seal  since  the  highest 
pressure  ratio  exists  at  that  pdnt,  independent  of  the  number  d*  knives  utilized. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

This  theas  has  presented  a  design  of  a  test  facility  to  aid  in  improving  sealing 
technology.  Based  upon  the  school's  available  facilities  and  the  levd  of  research  desired  to 
be  ccmducted  a  list  of  goals  was  presented  driving  the  design  of  the  test  rig.  The 
methodology  of  design  and  supporting  calculations  were  presented  alcmg  with  detailed 
mechanical  drawings  to  support  construction  of  such  a  test  rig.  A  proposed  data  acquiation 
system  was  presented  to  allow  adequate  collection  erf*  data,  and  initial  operational 
procedures  for  the  rig  were  included.  An  overview  of  the  factors  impacting  labyrinth  seal 
performance  were  presented  along  with  a  computer  model  to  predict  test  rig  operating 
parameters.  Graphs  of  the  test  rig  operating  realms  were  produced  along  with  a  required 
test  rig  parts  listing.  The  overall  conclusions  are: 

1 .  The  2D  test  rig  design  will  be  cost  effective,  simj^e  to  ccmstruct,  and  provide 
projected  performance. 

2.  The  proposed  design  met  the  design  goal  of  operating  with  a  factor  safety  greater 
than  3.3  at  conditions  of  Regap  <  100,000  and  Mach  <  0.9,  to  match  Reynolds  and 
Mach  numbers  that  occur  in  real  gas  turbines,  and  to  sup^rt  measurements 
significant  resolution  for  CFD  validation. 

4.  Modification  of  the  diffuser  assembly  may  be  necessary  due  to  the  lack  of 
information  concerning  dump-diffusers,  and  the  uncertainty  associated  with 
perf  oration  plate,  screen,  and  honeycennb  locations. 
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B.  RECOMMENDATIONS 

The  primary  recommendati(xi  based  upon  this  thesis  is  to  build  a  2D  test  rig 
according  to  the  provided  specifications.  Once  constructed  the  test  rig  should  be  used  to 
validate  the  MATLAB  program  flow  rate  and  chamber  pressure  predictions  by  experimental 
research,  and  eventually  modify  the  MATLAB  program  for  variable  carry-over  factors  and 
irregular  knife  gap  and  spacing.  The  rig  should  then  be  used  to  optimize  labyrinth  seal 
design.  After  the  basic  operation  of  the  rig  has  been  substantiated,  further  modifications  of 
the  rig  should  be  considered  to  improve  performance  as  follows; 

1.  Consider  an  eductor  system  to  create  a  vacuum  at  the  seal  exit,  this  would  provide 
additional  range  for  both  the  available  Reynolds  number  and  Mach  number  without 
raising  the  internal  operating  pressure  of  the  system. 

2.  Consider  perforations  and  suction  on  the  end  walls  to  minimize  boundary  layer 
growth  and  the  .subsequent  "end  wall  effect”. 

3.  Design  and  construct  a  3D  rotational  test  rig  to  determine  the  effects  rotation  on 
seal  performance  specifically,  this  rig  could  be  used  for  measurements  of  flow 
coefficients,  friction  factors,  and  heat  transfer. 

4.  Design  and  construct  a  modified  test  section  to  allow  testing  of  staggered  labyrinth 
seals. 
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STATOR 


Figure  3.1  Air  Supply  Facility  Diagram 


3.3  Manifold  Connection  Flange 


4.75"  FrqntView  9^25 


Figure  3.4  Diffuser  Assembly  Inlet  Plate 
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Figure  3.5  Diffuser  Assembly  Bottom  Plate 
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Figure  3.6  Diffuser  Assembly  Top  Plate 
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Front  View  ,  1"  overhang 

33.75"  \T 
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Figure  3.8  Diffuser  Assembly  Right  Side  Plate 
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Figure  3.9  Diffuser  Assembly  Perforation  Plates 


Side  View 

Note:  Both  plates  made  of  0.25" 

1st  stage  perforation  plate  Perforated  as  drawn.  2nd  stage  perforation  plate 
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Figure  3.10  Diffuser  Assembly  Perforation  Plates 
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Figure  3.1 1  Contraction  Section  Inlet  and  Exit  Flange 
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Settling  Section  Top&Bottom  plate  Top  View 
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Top  View 
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Top  View 
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Figure  4. 1  Data  Acquisition  Schematic 
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Figure  5.2  Flow  Rate  Correlation 


Figure  5.3  Chamber  Pressure  Correlation 
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Figure  5.4  Test  Rig  Pcrfonnancc 
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Figure  5.5  Test  Rig  Performance 
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Figure  5.6  Test  Rig  Performance 


Test  rig  performance,  Pexlt=35psia 
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Figure  5.7  TCvSt  Rig  Performance 
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Figure  5.8  Test  Rig  Performance 


jaquinu  spiouAsy 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

(O 

o 

m 

o 

u> 

o 

o 

CO 

CO 

CM 

CM 

T“ 

U) 

oinoinomoino 

^COCOCVICMt-;T^OO 

ooooooooo 


#MOM 


64 


Figure  5.9  Test  Rig  Performance 
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Figure  5. 10  Test  Rig  Performance 


3  knife  seal  45-15  psia 
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Figure  5. 12  Test  Rig  Performance 
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Figure  5.13  Te.st  Rig  Performance 
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Figure  5.14  Test  Rig  Performance 


APPENDIX  A:  SEAL  LEAKAGE  PREDICTION  PROGRAM 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  This  program  will  determine  the  1-0  intermediate  pressures  and  mass  flowthrough  % 
%  a  labyrinth  seal  with  2  knives.  The  user  specifies  seal  gap,  inlet  and  exit  pressures,  seal  % 
%  length  and  pressure  step  size.  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
clear 

%  Specify  seal  gap  in  inches. 
g=.i5; 

%  Specify  seal  length  in  inches  (Used  only  for  BETA  calculation). 

L=2: 

%  Specify  inlet  pressure  in  psia. 
po=75: 


%  Specify  outlet  pressure  in  psia. 
pn=i5: 

%  Specify  pressure  step  size  in  psia 
stepsize=2 

N=(po-pn)/stepsize 

%  This  converts  gap  and  length  to  meters. 
gm=g*.0254; 

Lm  =  L*.0254: 

%  This  converts  inlet  and  outlet  pressures  to  Pascals. 

pom=po*6894.76: 

pnm=pn*6894.76; 

stm  =  (pom-pnm)/N; 

%  This  specifies  gamma,  R  and  T  (in  Kelvin) 

R=287; 

T=293: 
mu  1  =.65; 


alpha  =  8. 52/(7. 23  +  (Lm/gm)): 
beta=i/((l  -alpha)''.5); 
mu=.65*beta; 


for i=1 ;N, 

pom(i,l)  =  pnm  +  (i*stm); 

p1(i,1)  =  sqrt(.5*pnm''2  +.5*((pom(i,l))''2)); 

q1  (i,1)  =  ((mul*gm)/sqrt(R*T))*sqrt((pom(i,l))"2-(p1  (i,1))''2); 
q2(i,l)  =  {(mul*gm)/sqrt(R*T))*sqrt((p1(i,l))''2-pnm''2); 
end 

pom 

pi 

pnm 

qi.q2 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  The  following  statements  modify  the  previous  program  for  seal  knives  of  3, 4, and  5.  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  This  is  for  3  knives. 

for  i=1 ;N, 

pom{i,1)spnm  +  (i*stm); 

pi  (i,1  )*sqrt((1/3)*pnm''2  +  (2/3)*(pom(i,  1  ))*2) ; 
p2(i,1  )=sqrt({2/3)*pnm''2  +  (1/3)*(pom(i,1))^2); 

q1(i.1)  =  ({mu1*gm)/sqrt(R*T))*sqrt((pom(i,1))>‘2-(pl  (i.1))*2): 
q2(i.1)  =  ((mul*gm)/sqrt(R*T))*sqrt({pl(i.i))^2-{p2(i.1))^2); 
q3(i,1  )  =  ((mu1  *gm)/sqrt(R*T))*sqrt({p2{i,1  ))''2-pnm'‘2); 
end 

%  This  is  for  4  knives, 
for  i=1 :N, 

pom(i,1)  =  pnm  +  (i*stm); 

p1(i.1)  =  sqrt(.25*pnm''2  +  .75*(pom(i.1  ))''2): 
p2(i,1  )=sqrt(.5*pnm'‘2  +  .5*(pom(i,l  ))'‘2); 
p3(i,  1  )»sqrt(.75*pnm'‘2  +  .25*(pom{i,1  ))''2): 

q1  (i,1  )  =  ((mu1*gm)/sqrt(R*T))*sqrt((pom(i.1))''2-(p1  (i,1  ))''2); 
q2(i.1)  =  ((mu1*gm)/sqrt(R*T))*sqrt((pl{i,1))^2-(p2(i,1))*2): 
q3(i,1  )*((mu1  *gm)/sqrt(R*T))*sqrt((p2(i.1  )K2-(p3(i,1))''2): 
q4(i,1  )  =  ((mu1  *gm)/sqrt(R*T))*sqrt((p3(i,1  ))^2-pnm''2); 
end 

%  This  is  for  5  knives, 
for  i=1 :N, 

pom(i,1)  =  pnm  +  (i*stm); 

pi  (i.1  )=sqrt(.2*pnm''2  +  .8*(pom{i,1  ))''2): 
p2(i,1  )=sqrt(.4*pnm*2  +  .6*(pom(i,1  ))''2); 
p3(i,1  )=sqrt(.6*pnm''2  +  .4*(pom(i,1  ))''2); 
p4(i,1  )=sqrt(.8*pnm''2  +  .2*(pom(i,1  ))''2); 

q1(i,1  )  =  ((mu1*gm)/sqrt(R*T))*sqrt((pom{i,1))''2-(p1  (i,1))''2): 
q2(i,1)  =  ((mu1*gm)/sqrt(R*T))*sqrt((p1  (i,1))''2-(p2(i,1))''2); 
q3(i,1)  =  ((mu1*gm)/sqrt(R‘T))*sqrt((p2(i,1))'‘2-(p3(i.1))''2): 
q4(i,1  )  =  ((mu1  *gm)/sqrt(R*T))*sqrt((p3(i,1))''2-(p4(i,1))''2): 
q5{i.1)  =  ((mu1*gm)/sqrt(R*T))*sqrt((p4(i.1))''2-pnm''2): 
end 
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Quantity 


20  FEET 


10  FEET 


Quantity 


APPENDIX  B:  LIST  OF  MATERIALS 


Airflow  piping «nd  valves 


Desciiptiao 


MANUAL  SHUT  OFF  VALVE  2*  GLOBE 


FULL  PORT 


ADJUSTABLE  PRESS  REG  VALVE 


(M5  PSIG,  FULL  PORT.  TIGHT  SHUT  OFF 


PRESSURE  RELIEF  VALVE,  ADJUSTABLE 


35-50  PSIG 


EXIT  VALVE.  4"  BUTTERFLY  OR 


GLOBE,  FULL  PORT 


FILTER  ASSEMBLY  GPI98-0030-020 


Diffuser  assembly 


DesciiptiOT 


1 "  STEEL  PLATE  9.25"X18.5" 


Approximate  cost 


(Dcdlais) 


$40 


Approximate  cost 


(Dollars) 


$30 


SCREEN/HONEYCOMB  ASSEMBLIES 
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Quantity 


Contraction  assembly 


Descripdan 


0.5-  STEEL  PLATE  9-X17- 


0.5"  STEEL  PLATE  18-X18- 


0.75"  STEEL  PLATE  1.5-X22" 


0.75"  STEEL  PLATE  2S"X22” 


1  -  STEEL  PLATE  0.75-X  10" 


0.375"  X3-  BOLTS 


Approximate  cost 


(DoUaus) 


$20 


Quantity 


2 


2 


2 


2 


18 


Setting  assembly 


Descripticn 


0.5"  STEEL  PLATE  2.5"X5" 


0.5"  STEEL  PLATE  5-X18" 


0.5"  STEEL  PLATE  5-X19- 


0.5"  STEEL  PLATE  4.5-Xl- 


0375"  X3"  BOLTS 


Ai^iroxiniate  cost 


(Dcdlars) 


$10 


Test  Section 


Quantity 


2 


1 


2  I 


_ i 

i 
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3 


Description 

Approximate  cost 

(Dollars) 

1.5"  STEEL  PLATE  2.5"X8" 

$30 

0.5"  STEEL  PLATE  8"X  18" 


0.5"  STEEL  PLATE  7"X  1.5" 


0.5"  STEEL  PLATE  1-X18" 


0.375"  X3"  BOLTS 


1-X4-X6-  glass  plate 
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APPENDIX  C:  FLOW  MODIFICATIONS 


SHEAR 
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